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ABSTRACT 

We assess the potential of nuclear starburst disks to obscure the Seyfert-like AGN that dominate the hard 
X-ray background at z ~ 1. Over 1200 starburst disk models, based on the theory developed by Thompson et 
al., are calculated for five input parameters: the black hole mass, the radial size of the starburst disk, the dust- 
to-gas ratio, the efficiency of angular momentum transport in the disk, and the gas fraction at the outer disk 
radius. We find that a large dust-to-gas ratio, a relatively small starburst disk, a significant gas mass fraction, 
and efficient angular momentum transport are all important to produce a starburst disk that can potentially 
obscure an AGN. The typical maximum star-formation rate in the disks is ~ 10 M yr _1 . Assuming no 
mass-loss due to outflows, the starburst disks feed gas onto the black hole at rates sufficient to produce hard 
X-ray luminosities of 10 43 -10 44 erg s _1 . The starburst disks themselves should be detectable at mid-infrared 
and radio wavelengths; at z = 0.8, the predicted fluxes are ~ 1 mJy at 24/im and ~ 10-30 fiJy at 1.4 GHz. 
Thus, we predict a large fraction of radio/X-ray matches in future deep radio surveys. Unfortunately, both the 
24/im and radio fluxes are comparable to those expected from the central AGN. In contrast, the starburst disks 
should be easily distinguished from AGN in future 100/im surveys by the Herschel Space Observatory with 
expected fluxes of ~ 5 mJy. Any AGN-obscuring starbursts will be associated with hot dust, independent of 
AGN heating, resulting in observable signatures for separating galactic and nuclear star-formation. This may 
be an explanation for the small observed L2-10 kev/ l 'L l ,(Q fiia) ratios observed from both z ~ and z ~ 1 
AGN. Finally, because of the competition between gas and star-formation, nuclear starbursts will be associated 
with lower-luminosity AGN. Thus, this phenomenon is a natural explanation for the observed decrease in the 
fraction of obscured AGN with luminosity. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — galaxies: Seyfert — galax- 
ies: starburst — X-rays: diffuse background 



1. INTRODUCTION 

All Active Galactic Nuclei (AGN) are powe red by accre- 
tion onto a central super massive black hole jLynden-Bell 
fl969HShakura & Sunvaevlll973HPringldfT98TtlBalbusll2003h . 
However, the observational manifestation of this accretion 
can vary dramatically from object to object. For example, 
the spectr al energy distrib utions (SED) of rapidly accreting 
quasars dElvis etal.ll 19941) are very different fr om the Galac- 
tic Center source Sgr A* (iMelia & Falckdl200ll) . or the black 
hole at the center of M87 (Ho 1999). In this case, the different 
observational properties are a result of the two distinct mech- 
anisms through which energy is liberated in radiatively effi- 
cient (for the quasars) and inefficient (for Sgr A* and M87) 
accretion flows. Another important parameter is the orien- 
tation of the black hole-accretion disk system to the line of 
sight. As a result, blazars, AGN that are viewed down the axis 
of a radio jet, show significantly di fferent SEDs a nd variabil- 
ity properties than all other AGN (iBach et al.ll2007l) . Finally, 
obscuration along the line-of-sight is another key parameter in 
determining the observational characteristics of an AGN. The 
vast majority of local opti cally-selected AGN show evidence 
for b eing obscured (e.g., iMaiolino & Riekd 119951: iHo et al.l 
119971) . and these are called Type 2 AGN, while the unobscured 
sources are defined as Type 1 AGN. Observationally, the ob- 
scuration is revea led as significant soft X-ray absorption (e.g., 
iTozzi et aill2006h. the disappearance of the o ptical broad per- 
mitted lines (iKhachikian & Weedmanl ll974). and a reddened 
continuum dWilkes et al.ll2.005b . 

Despite its near ubiquity, the origin of the obscura- 
tion around AGN is still unclear. Infrared e mission from 
dust near the sublimation temperature (e.g., iRiekd 119781: 



iNeugebauereTaTII 19791: lBarvainis|[T98H ISanders et al.ll 19891) 
and variability in the X-ray absorbing column density 
dRisaliti et al.l2002i 12005) both point to a location —1 pc from 
the central engine, although it is difficult t o rule out contribu- 
tions from the host galaxy at l arger scales (IMaiolino & Riekel 
119951: iBallantvne etai] 120031: iRigbv et al.ll2006l) . It is also 
possible that the origin of the obscuration may arise through 
a variety of processes for AGN with different luminosities o r 
at different points in their evolution (Ballanty ne et al.l l2006b). 

Perhaps the most significant challenge facing the pc-scale 
absorber (or obs curing 'torus' in the AGN unification model; 
lAntonuccill 19931) is to account for the observed 4-to-l Type 
2/Type 1 ratio of AGN r equired to fit the hard X-ray back- 
ground dGilli et al.l 12007). This fact requires the obscuring 
material to cover ^80% of the sky as observed from the cen- 
tral black hole, and thus have a structure so that its verti- 
cal scale h eight is of the sam e magnitude as its ra dius, i.e. 
H/r ~ 1. iKrolikl (120071) and IShi & Krolikl d2008h have ar- 
gued that the AGN emission absorbed and re-radiated in the 
infrared by dust in the torus can provide enough pressure sup- 
port to inflate a st ructure to the required scale. Similarly, 
Chang etaD, (120071) suggested that X-ray heating of the outer 
accretion disk may be sufficient to result in a Hj r ~ 1 geome- 
try. Another explanation, and the one considered here in more 
detail, is feedback from a nuclear starburst di sk (Fabia rTet al.1 
1 998 ; Wad a & Normanll200l iThompson et al.ll2005h . 

There are many reasons to consider starburst disks as a 
source of obscuration in AGN. At the most basic level, 
gas to fuel the black hole must be funneled toward the 
center of the galaxy and is therefore likely to cause star 
formation en route. In fact, recent adaptive optics ob- 
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servations by iDavies et al.l (120071) have shown the pres- 
ence of a post-starburst population at parsec scales in 9 
local Seyfert galaxies, indicating that this may indeed be 
a common phenomenon. In addition, recent models of 
AGN tori have suggested that the obscuring region may 
be clumpy (iNenkqya et alj|2002t [Dullemond & van Bemmel 
l200l iHonig et all 12006^ a c ommon trait of star-forming 
regions in starburst galaxies dForster Schreiber et al.l 120011 : 
Sniid ers et ail 120071) . Type 2 AGN also seem to be 
more correlated with star-formin g activity in general than 
the u n obscured Type 1 objects jC id Fernandes & Terlevich 
19951: lGonzaez Delgadoetaljj2001t iLevenson et al.l 120011: 
Kim et al1 l2006: Lacy et al. 2007). There is also evidence for 
a correlation between the observed intensity of star- formation 
in the host galaxy and the AGN nuclear luminosity dShi et al.l 
120071: IWatabe"lTaT]l2008l) . 

A more profound aspect of star-formation as a source 
of AGN obscuration is that it provides a connec- 
tion between the black hole- accretion disk envi r onmen t 
and the host galaxy (e.g., [Kawakatu & Wadal 120081) . 
It is now accepted that the correlations between the 
black h ole mass and bulge properties in early-type 
galaxies dM aggorian et a ll 119981: iFerrarese & Merrittl boOOt 
Geb hardt et al.ll2000t iTremaine et alj|2002ir are evidence for 
a significant connection between the grow th of the black hole 
and the build up of a galactic bulge (e.g.,[£ ilk & Ree s|ll998t 
Fabia nl 119991: iKaufmann & Haehneltl 120 00: Wv ithe & Loeb] 
2003[ iMurrav et alj|2005t iDi Matteo et alj|2005t iFabian et al.l 
2008). Therefore, it is expected that as the black hole and 
galaxy bulge are growing, there will be significant accretion 
on to the black h ole that is obscured b y star formation in the 
host galaxy (e.g.. iHopkins et all 1200 51). Indeed, the hard X- 
ray background, emitted by accreting black holes throughout 
the history of the Universe, has a very hard s pectrum, indica- 
tive of being dominated by obscured AGN dSetti & W oltier 
119891: iDe Luca & Molendil 12004 iGiih et al.ll2007l) . Interest- 
ingly, the redshift distribution of the obscured AGN that dom- 
inate the X-ray background peaks at z ~ 1, very similar to the 
peak in the cosmic star-formation history (|Tozzi et aT1l2001l: 
IBarger et al.ll200l lHopkinsll200lt IBarger et al.ll2005l) . This 
result may imply a connection between t he obscuring mate- 
rial and the evolution of the host galaxy ([ Fabian et al .1 119981: 
iFranceschini et al.l 119991: iBallantyne et all l2006al) . In fact, 
there is tentative evi dence that the fraction of obscured AGN 
does i ncrease with z (Ballanty ne et al.l2006atlTreister & Urrvl 
2006). Thus, if star-forming disks do provide a significant 
amount of AGN obscuration at z ~ 1 (where the X-ray back- 
ground sources are most common) then by studying the prop- 
erties and evolution of the absorbing material, we can directly 
probe the evolution of the underlying host galaxy. 

In this paper, we make use of the analytic starburst disk 
models developed by iThompson et al.l d2005l) to investigate 
the properties of starburst disks as a source of obscuring ma- 
terial in the Type 2 AG N found at z ~ 1. In the paper by 
Thompson et al.l (l2005h . the models were able to successfully 
describe many of the observed proprieties of Ultra-Luminous 
Infrared Galaxies (ULIRGs), and it was noted that under cer- 
tain conditions, the photosphere of the starburst disks may 
reach H/r ~ 1. Here, we adapt the model to consider 
the much less intense star-forming disks that might be ex- 
pected to be found around the Seyfert-like lumin osity AGN 
that dominate the popul ation at these redshifts (Ue da et al.1 
l2Q03t IBarger et al.ll2005l) . We search for the range of phys- 
ical parameters these starburst disks must have in order for 



them to provide significant obscuration. Both radio and far-IR 
fluxes are predicted to determine if the disks can be (or have 
already been) detected by current observatio nal surveys. The 
next s ection provides a brief review of the IThompson et al.l 
(120051) models, and then describes how the theory was al- 
tered and applied to this problem. Section |3]presents the re- 
sults of our calculations, and describes both the physical and 
observational properties of starburst disks that may obscure 
an average AGN. The results are discussed in Section 0] and 
we present our conclusions in Section When necessary, 
the following A-dominated cosmology is assumed in this pa- 
per: ff n = 70 km s" 1 Mpc -1 , fl A = 0.7, and il m = 0.3 
(ISpergel et al.ll2003l) . 

2. MODELS OF NUCLEAR STARBURST DISKS 
2.1. Review of Starburst Disk Models 

The properties of the nuclear starburst disks are calcu- 
lated using the one-dim ensional analytic theory developed by 
Thompson et al. (2005), and we solve the equations presented 
in their Appendix C. Here, we will briefly review the main as- 
sumptio ns and ideas behind the th eory, and refer readers to the 
paper by Thompson et al.1 (120051) for a complete discussion. 

The starburst disks are calculated as a single phase medium 
at a radius r from the center of a galaxy with a velocity disper- 
sion a and a supermassive black hole with mass Mbh- Fur- 
ther, the disk is assumed to be rotating at the Keplerian rate at 
all radii, and the gravitational potential is the sum of the point- 
mass of the black hole and an isothermal sphere describing the 
galactic bulge. The star-formation in the disk is a local phe- 
nomenon and is assumed to always adjust so that Toomre's 
Q parameter is unity; that is, Q = Knc s /7rGE g = 1, where 
kq is the epicycle frequency, fl is the Keplerian frequency, 
c s is the sound speed and S g is the gas surface density. This 
condition on Q allows the gas density p at any radius to be 
immediately written as a function of fi(r). The vertical sup- 
port of the disk is the infrared (IR) radiation pressure of the 
starburst against dust (the UV optical depth is always much 
greater than unity in these disks, so the vast majority of the ra- 
diation is re-radiated in the mid- and far-IR). Therefore, these 
starbursts are locally radiating at their Eddington limit. The 
vertical optical depth of the disk in the IR is r = £ g K/2, 
where k(T, p) is the Rosseland mean opacity. 

Fuel is fed to the outer radius of the disk at a rate of M ut- 
An unspecified global torque acts on the disk allowing the 
removal of angular momentum and gas to accrete slowly in- 
wards at a rate M. This torque is assumed to allow the radial 
velocity of the gas to be equal to a constant fraction m of 
the local sound speed; thus, M = 2irrY, g mc s . However, 
gas can be removed from the disk by local star-formation, 
so that the accretion rate through the disk can also be writ- 
ten as M — M ou t — J^° ut 27rr'S*c?r', where £* is the star- 
formation rate (SFR) per unit area. Therefore, there is a com- 
petition between star-formation using up the gas in the disk 
and the torque attemptin g to transport m aterial through the 
disk. As discussed in Thompson et alj d2005l) . this can be 
understood by considering the advection timescale through 
the disk, T ac i v = r/v r , and the star-formation timescale, 
t* = l/(?yf2), where v r is the radial velocity of the gas, and 77 
is the star-formation efficiency. If r ac i v < t* then material can 
move through the disk qu ickly enough that it can fuel a black 
hole while forming stars. Thompso n et alj 12005) showed that 
by equating these timescales at the outer radius i? out , one de- 
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rives a critical accretion rate, M cr it, below which nearly all 
the gas is converted into stars at large radii. 

The model disks are calculated from i? ou t inwards until a 
point is reached where dissipation by disk accretion domi- 
nates the heating over star-formation. At this point, Q = 1 
can no longer be maintained and the starburst disk structure 
calculation is terminated. The actual value of the inner ra- 
dius depends on the particular parameters of the model, but is 
typically < 1 pc. 

2.2. A Survey of Starbursts 

Thompso n et al.1 ((2005) calculated starburst disk structures 
appropriate for ULIRGs, but the X-ray background is dom- 
inated by AGN whose host galaxie s must have, on aver- 
age, S FRs much lower than ULIRGs (Ballantyne & Papovich 
2007). Thus, we have to calculate new starburst mod- 
els for this different scenario. Since it is clear that it is 
possible for starburst d isks to potentially obscure an AGN 
(iThompson et al.l 120051) . these new models are designed to 
cover a wide range of parameter space to determine the va- 
riety of properties that AGN-obscuring starburst disks may 
have. 

Each starburst disk calculation has 5 input parameters that 
were varied through different values: log(AfBH/Af©) = 7, 
7.5, 8, 8.5 (the galactic velocity dispersion is then calcu- 
lated through the Mr h-Q" relation: Mbh = 2 x 10 8 ct 4 ; e.g., 
iTremaine et al1l2002l) ; m = 0.0075, 0.01, 0.025, 0.05, 0.075, 
0.1, 0.2; R out = 50, 100, 150, 200, 250 pc; / gas = 0.1, 
0.5, 0.9 (th e gas fraction determin es M out and decreases to 
smaller r; Thomp son et al.1 |2005); and the dust-to-gas ra- 
tio of the gas (eit her lx, 5x or 10 x the local ISM value). 
As emphasized by Thom pson et al.l (120051) . the starburst disk 
structure is very sensitive to the temperature dependence of 
the Rosseland mean opacity. Similarly, the overall opac- 
ity value is an important parameter, as increasing the opac- 
ity can result in a larger temperature without necessarily en- 
hancing the SFR. There are numerous observational hints 
that the region around an AGN may have enhanc e d metallic- 
itv (iHamann & Ferlandlll999l: iNagao et all 120061; lArav et al 



l2007ir as well as increased grain growth (iChang et al.ll2007l 



and references therein). Therefore, it is important to calculate 
starburst disks with di ffering dust-to-gas ratios. Following 
Thompson et al.l (120051). we use the Ross eland mean opacity 
curve calculated by lSemenov et al.1 (120031) with the following 
parameters: a density of 10 6 cm~ 3 , "normal" silicates, homo- 
geneous distribution of dust materials, and spherical particles. 
To determine the opacity at different dust-to-gas ratios, this 
basic curve, appropriate for typical ISM value s, was multi- 
plied by the desired ratio (Ferguson et al. 20071). 

Starburst disk models were run for each permutation of 
the input parameters resulting in a total of 1260 distinct star- 
burst disk structures. However, before analyzing the prop- 
erties of these disk models, we must select the models that 
are most likely to both fuel and obscure a central active nu- 
cleus (i.e., ones where r a d v < t*). The potential for obscur- 
ing an AGN arises when the central temperature of the star- 
burst disk surpasses th e sublimation temperature of the dust 
(IThompson et al.ll2.Q05h . At that point, the opacity of the gas 
drops precipitously and the SFR rises dramatically in order to 
keep Q = 1. In this situation, the effective temperature of 
the disk (roughly the surface temperature) may remain below 
the sublimation temperature, leading to a sign i ficant vertical 
opacity gradient in the disk. IThompson et al.l (120051) argued 
that, although the disk may remain thin with a H/r < 1, the 
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FIG. 1. — Histogram of the radius where the SFR reaches its maximum. 
The models making up the histogram were selected so that the maximum 
star-formation rate occurred in the high temperature dust sublimation region 
that satisfied the two criteria described in the text. Of the 1260 starburst disk 
models calculated, a total of 512, or 41%, satisfy these conditions. The peak 
star-formation rate always occurs at r > 0.2 pc, with r Ri 1 and 3 pc being 
the two most common values. 

photosphere of the disk could puff up to H/r ~ 1. There- 
fore, each of the 1260 models was searched to determine if 
there existed a region that satisfied both of the following con- 
ditions: 

1 . the central temperature was greater than 900 K, 

2. the SFR was greater than 10% of the SFR at R out . 

The second criterion was necessary to exclude a small num- 
ber of models that marginally exceeded the temperature con- 
dition, but did not reach the sublimation temperature. Only 
512 models, or 41% of the calculations, produced an inner 
starburst that satisfied both of the imposed conditions. A 
histogram of the radius at which the SFR reached its maxi- 
mum within this region is shown in Figure Q] The histogram 
shows that the most common radius to find this burst of star- 
formation is r k 1-3 pc. Such radii are very similar to the 
scales of the obscuring medium mea sured in nearby AGN 
dJaffe et alj|2004t iTristram et al.ll2007l) . Thus, these starburst 
disk models reproduce the size scales necessary to obscure 
AGN. It is worth emphasizing that in choosing the range of 
model parameters, we made no a priori assumptions about 
what values would be more (or less) likely to produce these 
pc-scale starbursts, but only considered values to be typical of 
observed Seyfert-like galaxies. The distributions of observ- 
ables shown below are derived solely from the models shown 
in Fig. [TJ which we consider to be those that are most likely 
to obscure the central AGN. 

3. RESULTS 

3.1. Physical Properties of the Starburst Disks 

As a guide to the physical structure of one of the starburst 
disks, Figure [2] plots the radial dependence of several inter- 
esting variables from one representative model. The details 
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FIG. 2. — Plot of the mass accretion rate M (long-dashed line; in 
M© yr _1 ), the SFR M, (solid line; in Mq yr _1 ), the vertical optical depth 
t (dotted line), the central temperature T (short-dashed line; in K), and the 
effective temperature T a g (dot-dashed line; in K), all as a function of radius. 
The data are drawn from a starburst disk model with the following parame- 
ters: M B h = 10 s M Q> dust-to-gas factor= 1, m = 0.05, R ou t = 100 pc 
and / gas = 0.9. 

behind the models are described by Thomp son et al.l d2005l) . 
but it is worth emphasizing here that star-formation occurs 
throughout the disk (solid line), from i? ut = 100 pc to the 
inner-radius at < 1 pc. However, the central temperature of 
the disk (short-dashed line) increases to smaller radii as the 
rising density causes the disk to become more optically-thick 
to the stellar radiation. Eventually, this temperature reaches 
the sublimation temperature of the dust, and the opacity drops 
(dotted line), resulting in a significant burst of star-formation. 
Since the effective temperature (dot-dashed line) is still be- 
low this sublimation temperature, a significant vertical opac- 
ity gradient will exist in this region. Thus, this is the key 
region for where star-formation can puff up the disk photo- 
sphere to H/r ~ 1, and not at larger radii, even if the SFRs 
are comparable. This inner burst of star-formation nearly ex- 
hausts the gas supply in these models, but a small amount 
survives to make its way onto the inner accretion disk to feed 
the black hole (long-dashed line). 

We now turn to the physical properties for the ensemble of 
'successful' models shown in Fig.Q] What is necessary for 
these disks to be able to obscure an AGN on pc scales? The 
solid lines in Figure [3] show the distributions of the input pa- 
rameters derived from those models that have a pc-scale burst 
of star-formation. These plots show that successful models 
can be produced at every value of the input parameters, but 
pc-scale starbursts are more commonly generated if the disks 
have specific properties. For example, panel (a) shows that 
more massive black holes can more easily produce such star- 
bursts. This is perhaps unsurprising since in these Q = 1 
structures a larger black h ole mass results in a d enser and 
hence hotter starburst disk (Thom pson et al.ll2005l) . In con- 
trast, panel (b) shows that a smaller starburst disk is more 
likely to have the inner ring of star-formation. Smaller star- 



burst disks retain a larger fraction of their gas at a given radius 
than larger disks. Thus, the opacity is greater and the temper- 
ature can more easily reach the sublimation temperature of 
dust. The smallest i? ou t considered, 50 pc, is an interesting 
number as its roughly the same size as the star-forming re - 
gions found around nearby Seyferts by iDavies et alJ (2007). 
The models considered here do show that such small starburst 
disks are more likely to produce an AGN obscuring region 
than larger disks, in agreement with the observations. 

There is a similarly strong trend shown in panel (c) involv- 
ing the Mach number m. A larger value of m reduces r a d v , 
and thus causes more efficient mass accretion through the star- 
burst disk. Figure [5] shows that an inner burst of star forma- 
tion at pc-scales is more common at higher Mach numbers. 
In fact, ~ 50% of such models had m > 0.1. As with the 
smaller values of i? ou t, efficient mass accretion results in a 
larger fraction of gas at pc-scales, which increases the opac- 
ity and hence temperature of the gas. The same reasoning 
can also be applied to the dust-to-gas factor shown in panel 
(d), where clearly a larger dust-to-gas ratio results in a more 
common inner pc-scale burst of star-formation. A enhanced 
dust-to-gas ratio increases the vertical opacity to the starburst 
radiation, leading to larger temperatures and thus increasing 
the likelihood of surpassing the sublimation temperature. 

The gas fraction at i? out (Fig.Oe)) is another important pa- 
rameter for the existence of a pc-scale starburst, with f gas = 
0.1 producing only ~ 2% of the 'successful' models. The 
mass accretion rate at i? ou t is directly related to / gas , so that 
a larger accretion rate onto the outer disk is the result of in- 
creasing / gas . Hence, it is easier to exceed the Af cr i t crucial 
for feeding and obscuring the black hole. A larger / gas at i? out 
will also likely result in a larger fraction of gas on pc-scales, 
and the extra heating argument described above will hold. 
The accretion rate varies as f gas dThompson et alj|2005l) . thus 
models with / gas ~ 0.5 are adequate to produce potentially 
AGN-obscuring starburst disks. 

It is clear from the physics of the starburst model that in or- 
der to produce this pc-scale burst of star-formation that may 
obscure an AGN, the temperature of the disk must exceed 
the sublimation temperature of dust producing a sudden drop 
of the vertical opacity. Figure [3] shows that this situation is 
more likely to occur if one of the following is true: a high- 
mass black hole, a smaller starburst disk, efficient angular 
momentum transport in the disk, a reasonable gas mass frac- 
tion, or a large dust-to-gas ratio. However, not all of these 
conditions have to be simultaneously satisfied. For exam- 
ple, the dotted lines in Fig. [3] show the distribution of i? ut, 
/ gas and the dust-to-gas factor for the subset of models where 
log(MBH/-^©) = 7 and m = 0.025, two non-ideal con- 
ditions for the existence of the inner pc-scale starburst. In 
this case, the small starburst disks and large dust-to-gas ra- 
tios are basically required in order to generate a potentially 
AGN-obscuring starburst. Therefore, while it may be harder 
to produce a pc-scale burst of star-formation around a smaller 
mass black hole, it is not impossible if the other properties of 
the starburst disk are favorable. 

3.2. Star-formation Rates 

We now understand the conditions necessary for the star- 
burst disk to produce a potentially AGN obscuring ring of 
star-formation. Focusing in on the star-formation events 
themselves, Figure |4]plots the distribution of maximum SFRs 
found in these pc-scale starbursts. The plot shows that just 
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FIG. 3. — The different panels show the fraction of the models with pc-scale bursts of star-formation (i.e., those shown in Fig. [T} that have a specific input 
parameter. The panels show (a) the mass of the central black hole, (b) the outer radius of the starburst disk, (c) the logarithm of the Mach number m, (d) the 
dust-to-gas multiplicative factor and (e) the gas fraction at -Rout- The solid lines are the results when including all the successful models. The dotted lines plot 
the fractions when the additional constraints of logf-Meii/^©) = 7 and m = 0.025 are imposed. 
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FIG. 4. — Histogram of the maximum SFRs found in the models with 
pc-scale starbursts. About 55% of the models have a maximum SFR< 
20 Mq yr -1 . The distribution shows a weak tail to high SFRs, with a small 
number of models with peak SFRs > 300 Mq yr" 1 . The inset shows the 
same distribution, but now on a logarithmic axis. The most common SFRs 
are 10-30 Mq yr" 1 and only ~ 5% of the 'successful' calculations have 
SFRs > 300 Mq yr- 1 . 

over half the models have maximum SFRs < 20 M© yr -1 
(with values ~ 10 M yr -1 most common), r ates that are 
very typical of luminous IR galaxies (LIRGs; Wang et al. 
2006), but mu ch less than those observed in ULIRGs (> 
100 M Q yr -1 ; iRodriguez Zaurfn et all 120081) . However, the 
star-formation per unit area in the model bursts are very 
high (~ 10 6 M Q yr -1 kpc" 2 ), excee ding those observed 
in UL IRGs (~ 10 3 M Q yr^ 1 kpc" 2 ; iGracia-Carpio et all 
2008). Thus, these pc-scale starbursts are basically minia- 
ture ULIRGs: small regions of very efficient star-formation, 
but because the area is so s mall, the maximum SFRs can 
be LIRG-like. Interestingly, Bal lantyne & Papovichl (120071) 
showed that the host galaxies of the AGN that dominate the X- 
ray background must, on average, have SFRs< 100 M Q yr -1 
to avoid violating the cosmic IR background. This constraint 
is consistent with the predictions of Fig. @] where ~ 80% of 
the models have maximum SFRs less than 100 M Q yr -1 . 

The distribution of SFRs shows a small number of mod- 
els with maximum SFRs > 300 M Q yr -1 . The Fig. |4] inset 
indicates that only 5% of the models present such ULIRG- 
like SFRs. These massively star-forming disks only occur 
when log(M/M Q ) = 8.5, m = 0.1 or 0.2, and / gas > 0.90. 
Clearly, significant gas resources are necessary to power such 
bursts of star-formation. 

3.3. Observational Consequences 
3.3.1. AGN Properties 



As mentioned in Sect. 13. II if T a d v < r* a small amount 
of gas survives its passage through the starburst disk and will 
feed onto the sub-pc accretion disk, and may ultimately ac- 
crete onto the black hole. Under the assumption that the 
accretion rate at the inner-most edge of the starburst disk is 
the same as the accretion onto the black hole, the bolomet- 
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FIG. 5. — Histogram of the Eddington ratio, or ibol/^Edd. predicted 
from a central AGN being fed from the surrounding starburst. Mass-loss in 
the inner accretion disk is not included, so these ratios are strict upper-limits. 
Sixty percent of the models have 0.1 < ibol/^Edd < 1. with a further 
~ 20% apparently accreting at super-Eddington rates. 



ric luminosity of the resulting AGN can be estimated from 
Lbo\ = r)Mc 2 , where r\ = 0.1 is the radiative efficiency of 
the accretion disk. This is not a very realistic assumption, 
since outflows from the central regions of AGN are commonly 
observed, as well as being a standard pr ediction from numeri- 
cal simulations of acc retion flows (e.g. JCrenshaw et "ai]|2003t 
lHawlev & K rolik 2006). However, these outflows ensure that 
the accretion rate onto the black hole is always less than the 
accretion rate at the outer edge of the accretion disk, thus the 
luminosities calculated from the predicted M are strict upper- 
limits 1 . A further, implicit, assumption in this calculation is 
that there are no other sources of fuel for the black hole aside 
from that provided by the starburst disk. 

Keeping those assumptions in mind, we plot in Figure [5] 
the Eddington ratio Lboi/^Edd distribution from the models 
with pc-scale starbursts. Here, L^dd = &KGMwm v c/ o~t 
is the Eddington ratio for spherical accretion, where m p is 
the proton mass and <tt is the Thomson cross-section. The 
histogram shows a log-normal-like distribution with a peak 
at iboi/^Edd ~ 0.5. In all, about 20% of the models have 
Eddington ratios less than 0.1, ~ 60% have ratios 0.1 < 
iboi/iEdd < 1, while ~ 20% predict super-Eddington accre- 
tion. Since these ratios will be upper-limits, we find that the 
majority of the Eddington ratios predicted by fueling from the 
model starburst disks are consis tent with the valu es inferred 
from obscured, z ~ 1 AGN (e.g. lBallo etai1l2007l) . 

While the Eddington ratio is an interesting parameter, it 
is not a true observable. Assuming a standard radiatively- 
efficient AGN spectrum, bolometric corrections can be used 

1 The mass-loss rate from an accretion disk depends on the location from 
where it is launched and how it is driven, bo th extremel y uncertain numbers 
for almost all AGN ICrenshaw efafl feOfB). Recently, Mil ler et alj f200Sl) 
argued that an outflow from the Galactic black hole GRO J1655-40 was 
launched from the disk by magnetic processes. The mass-loss rate in this 
case was estimated to be > 5-10% of the accretion rate. 
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FIG. 6. — Histogram of the 2-10 keV luminosity of AGN fed from the 
starbursts disks with pc-scale bursts of star-formation. These are strict upper- 
limits. The range of luminosities is consiste nt with those AGN observed to 
dominate the X-ray background at z ~ 1 lUeda et al. 2003; Barger et al. 
12001) . 

to determine th e AGN luminosity in the 2 -10 keV X-ray 
band. Recently, IVasudevan & Fabian] (|2007) determined the 
2-10 keV bolometric correction for a number of nearby AGN 
with both X-ray and UV data. These authors found that the 
bolometric correction may depend on the Eddington ratio of 
the AGN, with a value of L^ D \/Lx ~ 20 for Lboi/^Edd < 
0.1, while L^oi/Lx ~ 55 for larger values, where Lx is 
the AGN luminosity in the 2-10 keV band. Employing these 
bolometric corrections, strict upper-limits to Lx can be cal- 
culated for our models and its distribution is plotted in Fig- 
ure [6] Over 90% of the predicted X-ray luminosities fall be- 
tween Lx = 10 43 erg s _1 and 10 44 5 erg s~ 4 , which is pre- 
cisely the range of luminosities observed from th e AGN that 
dominate the hard X-ray background at z ~ 1 ( Ue da et all 
2003; Bar ger et af] 12005). It is interesting that although the 
nuclear accretion process is not modeled, the AGN accretion 
rates predicted by the competition between star-formation and 
a global feeding process on scales of tens of pc gives such 
good agreement with the observations. Although the details 
are unclear, this result may support models of AGN feeding 
through s tar-f orming disks and global, bar-like instabilities 
(see Sect. g3]i. 

3.3.2. The Starburst Disk 

To determine the observational properties of the starburst 
disk itself, we calculate a simple SED for each model assum- 
ing bl ackbody emission at ea ch radius and a face-on viewing 
angle dThompson et al. 2005): 



XL X = 



2irhc 2 



lixrdr 



e^qy[hc/Xk B T cS (r)} - 1 ' 



(1) 



where R[ n is the inner radius for each model, h is Planck's 
constant and fee is Boltzmann's constant. An example SED, 
calculated from the same starburst model as described in 
Fig. |2] is shown in Figure [7] The SED shows two maxima: 



M BH =10« M 

dust/gas factor = 1 

m=0.05 

R out =100 pc 
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FIG. 7. — Example SED of a starburst disk from a model with the fol- 
lowing parameters: A/bh = 10 s Mq, dust-to-gas factor= 1, m = 0.05, 
Rout = 100 pc and / gas = 0.9. The inner pc-scale starburst causes a bump 
in emission at A ~ 5/xm. The SED is calculated from eq. [T]and assumes a 
face-on geometry; hence the luminosity is an upper-limit. 



one at ~ 30/im from the star-formation at the outer edge of 
the disk, and one at ~ 5/im from the pc-scale starburst. This 
shorter-wavelength peak is a generic feature of the models 
with the inner ring of star-formation, and is due to the loss of 
dust opacity in this region that limits the reprocessing avail- 
able for the hot starburst emission. Thus, the effective temper- 
ature at these radii moves into the near-IR regime. As these 
SEDs are calculated assuming the disk is viewed face-on, the 
luminosities and fluxes discussed below are upper-limits. 

Figure[8]plots the distribution of 24/im flux at z — 0.8 pre- 
dicted from the starburst disk SEDs. The choice of z = 0.8 
is illustrative, but is a common redshift for the population of 
obscu red AGN that we wish to describe (e.g., iBarger et alj 
2005). The most common 24 jim flux predicted by the star- 
burst models is ~ 1 mJy, but values can be as low as 10 /iJy, or 
as high as 30 mJy. As mentioned above, these values should 
be treated as upper-limits. The blue shaded region shows 
where the models with a maximum SFR< 30 M Q yr~ 4 , 
while the red area plots models with a maximum SFR> 
100 M yr . The starburst disks with lower SFRs, make 
up over 50% of the models (Fig. [4j and dominate the 24 /jm 
predictions at fluxes < 1 mJy. For fluxes between 1-3 mJy, 
the models with intermediate SFRs dominate before ceding to 
the small number of high SFR disks. 

The fluxes shown in Fig. [8] are well within the sensitivity 
limits of many of the deep Spitzer surveys over the last several 
years dLacv et al . 2004; Rigby et al. 200 41; iFranceschini et al J 
2001 IStern et al l 120051: lAlonso-Herrero et alj 12006 



Barmby et al. 2006; Brand et al. 2006; Perez-Gonzalez et al. 



2008). Thus, if these disks exist around AGN, they should 
have been detected by these surveys. However, the AGN 
itself can heat dust to near- and mid-IR emitting temperatures. 
The calculations presented in this paper do not include the 
feedback of the central AGN on the starburst disk. To 
estimate the 24 /im flux from AGN heating, we make use of 
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FIG. 8. — The solid line plots the distribution of predicted 24/im flux at z = 
0.8 from models with pc-scale starbursts. The blue region denotes the loca- 
tion of models with a maximum SFR < 30 M yr" 1 , while the red region 
shows the distribution of models with a maximum SFR > 100 Mr yr — . 
All of these fluxes are derived from SEDs calculated by equation [TJand are 
thus upper-limits. The dashed and dotted lines are predictions of the 24 (im 
flux due to AGN heating of a constant den sity torus at a distance of 1 and 
10 pc, respectively I Ballantyne et al. 2006b). These calculations made use of 
the distribution of X-ray luminosities from Fig. [6] 
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FIG. 9. — The solid histogram plots the distribution of 1 .4 GHz radi o fluxes 
at z = 0.8 as estimated from the radio-IR correlation of Bell (2003). As in 
Fig. [8] the blue and red regions denote the location of the models with max- 
imum SFRs < 30 and > 100 Mq yr -1 , respectively. As these radio fluxes 
were derived from the starburst SEDs, they should also be treated as upper- 
limits. The dashed histogram is an estimate of the nuclear radio-emission 
from the central AGN assuming it is radi o quiet, has a flat spectrum of a = 
and vL v {b GHz) = 10 _5 L 2 -io keV fTerashima & W ilson 2003). 



the calculations by Ball antyne et al.l d2006bl) who modeled 
AGN SEDs from the X-ray to the far-IR assuming a constant 
density torus with no star-formation. These models allow one 
to easily connect the X-ray luminosity to predicted IR fluxes, 
so the distribution of X-ray luminosities shown in Fig. [6] can 
then be converted to a distribution of 24/im flux with the 
results plotted in Fig. [8] The dashed and dotted lines assume 
a distance from the AGN to the inner radius of the torus 
of 1 and 10 pc, respectively. The lines were derived from 
models that are averaged over all X-ray column densities, and 
assume a (1 + z) 03 evolu t ion of the AGN Type 2/ Type 1 ratio 
dBallantyne et al] l2006ah iTreister & Urryl |2006), although 
this assumption does not affect the results shown here. We 
see that a dust torus heated by a distribution of AGN X-ray 
luminosities from Fig.|6]results in very similar 24/im fluxes as 
that predicted from the starburst disks. Thus, we recover the 
long-standing problem of attempting to determine the relative 
influence of star-formation and AGN heating in an observed 
spectrum. However, the nuclear starburst models described 
here present a new wrinkle: observable hot-dust emission in 
the near- and mid-IR that is not caused by AGN heating. This 
may mislead some methods of AGN identification in the IR 
that a re based on the shap e of the near- and mid-IR spectrum 
(e.g., iDonlev et al.l |2007). We conclude that the nuclear 
starburst disks will be easily visible in the mid-IR, but since 
the predicted fluxes are nearly identical to those expected 
from AGN heating, it will be extremely difficult to positively 
identify them or measure their properties (see Sect. |4)>. 

Another observational consequence of the nuclear starburst 
disks will be radio emission. There have been numerous at- 
tempts to link th e properties of AGN found in deep radio and 
X-ray surveys (Bau er et alJ l2002t iGiorgakakis etal] l2004t 



Simpson et al. 2006; Barger et al. 2007; [Richards et d]|2007t 
Rovilo s et al. 2007|; lPark et al.l2008l) . but the matched sources 
seems to only c omprise a small fraction of all AGN (e.g., 
Barger et al. 2007). Radio fluxes at 1.4 GHz can be predicted 
for the starburst disk by m aking use of the IR SEDs and the 
radio-IR correlation 2 (e.g.. lYun et al.ll200ll) . We spec i fically 
use the form of the radio-IR correlation derived bv lBelll (120031) 
that employs the total 8-1000 /im luminosity to make sure we 
account for the emission f r om th e pc- scale starburst. Obser- 
vations by App leton etail d2004l) and llbar etal] (120081) have 
shown that the l ocal radio-IR correl ation does seem to hold at 
high z (but see lBeswick et al.l l2008). Therefore, up on assum- 
ing a typical starburst radio spectrum of a = 0.8 (lYun et alJ 
1200 U where S v oc v~ a ), we plot in Figure [9] the distribu- 
tion of radio fluxes from the 'successful' starburst disks at a 
representative redshift of z = 0.8. As in the previous figure 
showing the 24/im predictions, the contribution from mod- 
els with a maximum SFR < 30 and > 100 M Q yr" 1 are 
shown by the blue and red hatched regions, respectively. The 
radio flux estimates are again upper-limits as they originate 
from the starburst IR SEDs. The most common radio flux 
predicted for these starburst disks is ~ 10-30 /jjy, values 
that will again be dominated by disks with maximum SFRs 
< 30 M Q yr _1 . This is an interesting flux because it is 
just below the sensitivit y limit (~ 40-50 /j,Jy) of the deep - 
est radi o surveys (e. g., Bauer et al. 20021 iBarger e t al. 2007). 
Indeed, iBauer et ail d2002l) found that the fraction of X-ray 
AGN with radio matches increased as they approached their 

2 The assumption that such intense starbursts will follow the observed 
radio-IR correlations may be problematic, as this might require very large 
magnetic energy densities (e.g., Thompson et al. 2006). This problem will be 
investigated in future work. 
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sensitivity limit, and suggested that deeper observations will 
find a higher fraction of matches. This conclusion agrees with 
the predictions of Fig. [9] 

Inte restingly, fGeorgaka kis et al.l d2004l) and iRichards et alj 
(120071) find a slight tendency for the fraction of radio/X-ray 
matches to increase with X-ray obscuration. That is, the more 
obscured an AGN, the more likely that it will have a radio 
counterpart i n the deep radio sur veys. This result, although 
still tentative (Rovilos et al. 2007), is consistent with our pro- 
posed scenario where obscured AGN should be associated 
with significant star-formation. Future work is needed in the 
modeling to connect the starburst disks to predictions of X-ray 
column densities. Further deep radio observations are also re- 
quired to confirm this trend. 

As with the 24/im flux, a major observational difficulty 
is determining whether the origin of the radio flux is due 
to the AGN or to star-formation. It has long been known 
that the radio source population seems to change at a flux of 
~ 1 mJy, with radio-loud AGN dominating the counts above 
this flux, and, it was originally thought, star-forming galaxies 
dominating below this flux (ICondonlll984t I Windhorst et alj 
19851; ISevmour et alj 120041) . However, as suggested by 
Jarvis & Rawlin gs (2004), radio quiet AGN will start to con- 



tribute at faint radio fluxes, and, indeed, recent observations 
find that AGN contribut e ~ 50% of the radio population 
at the ~ 50 iiJy level (|Simpson et al.l 120061 : ISevmour et alJ 
120081: ISmolcic et all 120081) . To estimate the radio-flux due 
to the nuclear emission from the AGN, we use the result 
of iTer ashima & Wilsonl (120031) who showed that for local 
radio-quiet Seyfert galaxies and quasars, vL v (h GHz) = 
10~ 5 Zv2-io kcV- Combining this result with the upper-limits 
for the X-ray luminosities plotted in Fig. [6] results in a distri- 
bution of radio flux due to the AGN itself (dashed histogram 
in Fig.[9]l. To calculate the 1.4 GHz flux, a flat a — spec- 
trum was assumed, typical of the cor e nuclear emission from 
an AGN dKellermann & OwenlfT988l) . We find that the AGN 
radio flux at z — 0.8 are predicted to be only a factor of 
~ 3 lower than the starburst flux, indicating that confusion 
between AGN and star-formation will be a factor in studying 
this population. However, if the data are available, methods 
such as spectral indices, radio morphology and radio-to-mid- 
IR flux ratios have been shown to succe ssfully discriminate 
between the two po ssible origins (e.g.. iDonley et all 2005: 
ISevmour et aHl2008l) . 

As noted above, both 24/im and 1.4 GHz observations of 
distant AGN will have difficulty distinguishing the presence 
of a nuclear star-forming disk from the AGN emission. The 
starburst SEDs show a wide range of emission because the 
star-formation occurs over a broad range of temperatures. 
Therefore, it is useful to predict how these starbursts disks 
may appear in the far-IR, especially with the imminent launch 
of the Herschel Space Observatory. Figure QJJ plots the 
100/im and 350/im flux distributions from our starburst mod- 
els at z — 0.8. Again, since these predictions are derived 
from the SED calculated using eq. [JJ they should be consid- 
ered as upper-limits. Both of t hese distributions predict fluxes 
~ 1-10 mJy. The models of Ball antvne et al.l j2006b) were 
employed to estimate the contribution from an AGN heated 
torus at these wavelengths. At both these wavelengths, the 
torus emission was predicted to be over an order of magnitude 
smaller. Those calculations assumed a relatively compact ge- 
ometry for the AGN absorber, it is possible a more extended 
Compton-thick absorber could produce greater far-IR emis- 
sion, but we defer a calculation of such a scenario until later 
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FIG. 10. — The solid histogram plots the distribution of 100/tm flux pre- 
dicted from the starburst disk models at z = 0.8. The dashed line plots the 
same results at 350 /im. 

work. 

Two Herschel key projects are the PACS Evolutionary 
Probe 3 that has a planned sensitivity limit at 100/im of 1- 
6 mJy. At 350^m the Herschel Multi-tiered Extragalactic Sur- 
vey 4 (HerMES) predicts a sensitivity of 10-20 mJy. Thus, we 
predict that the PACS Evolutionary Probe survey should de- 
tect these starburst disks associated with hard X-ray sources. 
The HerMES survey will likely only detect the very rare high 
star-forming or low-redshift sources. 

To summarize, the nuclear starburst disks proposed in this 
paper may reveal themselves to observers in multiple ways. 
They should have already been detected in the mid-IR by 
Spitzer surveys, but the predicted fluxes are of the same or- 
der as those from AGN heating of a non-star-forming torus, 
leading to the familiar IR spectral decomposition problem. 
At 1.4 GHz, the starburst disks are predicted to be slightly 
fainter than the deepest current radio surveys. Thus, a s al- 
ready indicated from previous work dBauer e t al. 2002), the 
fraction of radio/X-ray matches should increase as surveys 
probe to fainter flux levels. Unfortunately, the predicted radio 
flux levels are comparable to what is expected from the cen- 
tral radio-quiet AGN. Determining the radio spectral index or 
radio-to-IR flux ratios are two promising techniques that have 
been shown to distinguish bet ween star-formation a nd nuclear 
activity in radio sources (e.g.. ISevmour et aLll2008l) . Finally, 
the cleanest test may occur in the far-IR, where the contami- 
nation from AGN heating is expected to be minimized. Fig- 
ure [lO] shows that these disks should be easily detectable by 
the planned Herschel surveys at 100 /im. The question then 
arises: if star-formation is detected in these AGN host galax- 
ies, is there a way to determine if it is originating in the nu- 
cleus? We address this question in the next section. 

4. DISCUSSION 

3 http ://w w w. mpe.mpg . de/ir/Research/PEP/ 

4 http://astronomy.sussex.ac.uk/~sjo/Hermes/ 
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4. 1 . Distinguishing Between Nuclear and Galactic 

Starbursts 

In the previous section, some of the observational prop- 
erties of the nuclear starburst disks were explored, and we 
showed that they should be detectable at both radio and IR 
wavelengths. However, at the redshifts of interest (z ~ 1), 
it is more than likely that star-formation will be ongoing at 
galactic scales within the AGN host galaxy. Thus, unless the 
galaxy is resolved on a kpc scale or less, any observational de- 
tection of star-formation signatures from an AGN host galaxy 
may be ambiguous with regard to its location. 

Yet, there is a key difference between the nuclear starbursts 
presented here and the more usual star-formation that is on- 
going on large scales. As described in Sections |2. 1 1 and [37X1 
the crucial pc-scale burst of star-formation only occurs when 
the central temperature of the disk exceeds the dust sublima- 
tion temperature. This region of the disk does not reach these 
temperatures by heating by an AGN, but rather as result of 
the increasing density required by a star-forming Q = 1 disk 
in a black hole potential. Thus, in contrast to typical Galac- 
tic star-forming regions, these models predict significant star- 
formation from a region with very hot dust emission, but not 
AGN heated dust. Therefore, if such nuclear starbursts are 
common, there should be more near-IR or 24/im flux than 
expected from just pure AGN heating. As mentioned above, 
if this is true then it may complicate AG N selection based o n 
power-law like or hot mid-IR colors (e.g jDonlev et al.l2007l) . 

There may be evidence for this effect already in the lit- 
erature. Several observations of both nearby and high- 
redshift AGN have shown that the ratio of the absorption- 
corrected Lx to a near- or mid-I R luminosity is relatively con- 
stant with column densi ty Nn dLutz et al.ll2004t iRigby et all 
l2006t iHorst et aT]|2008l) . This is not expected from simple 
AGN unification and torus models where the observed mid- 
IR luminosity is pr edicted to depend significantly on Nn 
dPier & Krolikl fT993). Evidently, there is still hot dust emis- 
sion observed from AGN at all column densities. This fact 
has been used to argue that the AGN torus should be clumpy 
(IHorst e t al. 2008), but it is also consistent with a nuclear star- 
burst disk around the AGN. To show this we plot in Figure [TTI 
the distribution of the £2-10 kcV /vL v (Q /im) ratio for the 
starburst models. This histogram should only be taken as in- 
dicative, since both the X-ray and IR luminosities are upper- 
limits. However, taking the results at face-value, we find that 
the peak o f the distribution is ve ry similar to the va lues ob- 
served by iLutz et al.1 (120041) and IRigbv et all (120061) Inter- 
estingly, the AGN torus models of Ballanty ne et al.l |2006b) 
were not able to reproduce these ratios, as they we re not able 
to pro duce a sufficient near-IR luminosity. While ILutz et all 
(2004) argued that they corrected their data for star-formation 
based on PAH emission, if the starburst is produced in these 
hot pc-scale reg ions, the PAH emission may be overwhelmed 
dDesai et al.l2007l) and the star-formation rate underestimated. 
Thus, nuclear star-formation may be another mechanism to 
explain the ubiquitous presence of hot dust emission even at 
large column densities. 

4.2. Obscuration as a Function of AGN Luminosity 

Another interesting aspect of AGN phenomenology is the 
observed decrease in the f raction of obscured AGN with in- 
creasing AGN luminosity (lUeda et alJl2003HLa Franca et al.1 
120051; lAkvlas et al.ll2006t iTreister et al.ll2008l) . This observa- 
tion has been most commonly explained by the receding torus 



0.3 



0.2 



0.1 



log(L 2 





v /vL (6/mi)) 



FIG. 11. — Histogram of the estimated L2— 10 keV /^L^ft /im) ratio 
predicted for the nuclear starburst models. They peak i n the exact same 
range as the obser vations from both local and 2 ~ 1 AGN iLutz et alJ2004l : 
IRigbv et aU2006l) . 



model dLawrencd 1 99 U [Simpson 2005) which argues that the 
expansion of the dust sublimation radius with luminosity will 
push back the torus to larger radii and therefore decrease its 
covering factor. In this paper, we propose that the AGN ob- 
scuration is dominated by a nuclear starburst at a distance of 
~ 1 pc. In this section, we investigate if this model can also 
provide an explanation for the variation of the Type 2/Type 1 
ratio with luminosity. 

Figure [12] plots how two of the starburst parameters de- 
pend on the estimated X-ray luminosity. The left panel of 
the figure plots the distribution of the maximum SFRs found 
in the model starburst disks. The blue and red hatched 
regions denote the contributions from models with 42 < 
log(Lx/erg s -1 ) < 43.5 and 44 < log(L x /erg s _1 ) < 
44.5, respectively. This plot shows that the most luminous 
AGN are associated with the most intense star-f ormation, 
consi s tent with recent ob servational suggestions (IShi et alJ 
120071; IWatabe et alj|2008l) . In contrast, the lowest luminos- 
ity AGN are predicted with starburst disks with only mod- 
est star-formation rates. The right panel of Figure [121 plots 
the distribution of the input angular momentum parameter m. 
In this case, the dotted and dashed line denote the contribu- 
tions from models with 42 < \og(Lx /erg s _1 ) < 43.5 and 
44 < log(Ljs:/erg s _1 ) < 44.5, respectively. As with the 
maximum SFRs, these two ranges of X-ray luminosities are 
separated in the diagram with large values of m only associ- 
ated with luminous AGN. Similarly, starburst disks with weak 
angular momentum transport (i.e., low values of m) seem to 
only be able to feed a low-luminosity AGN. 

These results can be used to posit an explanation for the ob- 
served decrease in AGN obscuration fraction with luminos- 
ity. Recall that in this starburst disk model there is a com- 
petition between star-formation and gas accretion through the 
disk. If there is intense star-formation within the disk then 
this reduces the gas supply available to fuel a black hole. A 
circumnuclear starburst will therefore more likely be associ- 
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FIG. 12. — (Left) The solid histogram plots the distribution of maximum 
SFRs seen in the models (as seen in Fig. [4}- The blue and and red hatched re- 
gions denote the contributions from models with 42 < log(Lx/erg s —1 ) < 
43.5 and 44 < log(Lx/erg s™ 1 ) < 44.5, respectively. (Right) The 
solid line plots the distribution of the angular momentum parameter m 
found in the 'successful' starburst disk models (as seen in Fig. [3jc)). The 
dotted and dashed lines denote the contributions from models with 42 < 
log(Lx/erg s _1 ) < 43.5 and 44 < log(Lx/erg s _1 ) < 44.5, respec- 
tively. 

ated with both a lower-luminosity and obscured AGN. As the 
luminosity of the AGN increases, it is necessary to also in- 
crease the gas supply onto the black hole. This could be done, 
for example, by more efficient angular momentum transport 
in the starburst disk (Fig. IT2Ti. While higher luminosity AGN 
can be accompanied by very intense nuclear starbursts that 
may also obscure the AGN, such bursts will likely be short- 
lived. More luminous AGN, such as quasars, are unlikely 
to be fueled through these starburst disks, because circumnu- 
clear star-formation would consume the necessary gas supply. 
The obscuration responsible for Type 2 quasars would then be 
generated at larger, galactic scales a s a result of the pow erful 
events shaping the host galaxy (e.g. jHopkins et al. 2005). 

This starburst scenario has an advantage over alternative 
mechanisms that r ely on AGN radiation to produce an obscur- 
ing medium (e.g., lKrolikll2007t IChang et al.1l2007l) . In these 
situations, more obscuration is predicted with a larger AGN 
luminosity, as the vertical pressure support of the obscuring 
medium is directly related to the intercepted AGN luminosity. 
A significant outflow will need to develop in order to limit the 
covering factor. While it is clear that a more detailed model of 
nuclear starburst disks is required to make quantitative predic- 
tions, we find that, at least qualitatively, identifying a pc-scale 
starburst as the origin of the obscuration for some AGN can 
provide a natural explanation for the observed decrease in the 
fraction of Type 2 AGN with luminosity. 

4.3. Torques and Angular Momentum Transport 

One of the most important parameters in determining 
whether a circumnuclear starburst can both feed and obscure 
an active nucleus is the efficiency of angular momentum trans- 
port through the disk. In this paper we have assumed that a 
global torque, such as a bar or spiral waves, operates on the 
starburst and removes angular momentum. We further as- 
sume that the resulting accretion rate through the disk can 
be sustained at a Mach number m. The results of Sect. [3] 



show that the likelihood for the formation of a pc-scale star- 
burst that could obscure an AGN is very sensitive to m, with 
values of 0.1-0.2 being most favorable, but lower values 
are certainly pos sible, especially for lower-luminosity AGN. 
iGoodmanl (2003 ) argues that spiral waves in an accretion disk 
could cause inflow with m ~ 0.1, but it is not clear whether or 
not such global coherent structures will be common for these 
disks. Some insight may b e gleaned fr om the study of gas 
flows in galactic nuclei (see Wada 2004 for a recent review). 
High-resolution simulations of the gas dynamics in a nuclear 
potential are now regularly finding coherent structures such 
as ba r s or spirals in the gas flows dEnglmaier & Shl osman 
2004; Maciejewski 2004). These structures do remove an- 
gular momentum and funnel material to the galactic cen- 
ter, although the Mach number is not specified. These gas 
flows may ultima t ely be come turbulent which, as argued by 
Wada & Normanl §002), may also cause further gas accre- 
tion. 

Alternatively, a local visco sity mechanism, such a s the 
magneto-rotational instability (Balbus & Hawley 1998), may 
operate in the starburst disk. To determine how this differ- 
ent angular transport process might affect the development of 
a pc-scale starburst, calculations of star burst disks were run 
with a simple a viscosity prescription (IShakura & Sunvaevl 
119731: see also Thomp son et al.l 120051) . For simplicity, a 
value of a — 0.3 was considered with all the other param- 
eters (Mbh> ^out> /gas, and the dust-to-gas factor) varying 
as before. Fifty percent of the 180 models generated a pc- 
scale starbursts with the parameters exhibiting very similar 
trends to the on e s show n in Fig. [3] However, as found by 
iThompson et al.l (120051) . this local viscosity does not move 
material inward fast enough, and the resulting mass accretion 
rate remaining to fuel the black hole is much smaller in this 
scenario than in the global torque models described above. In 
fact, over half of the 'successful' a-disk models resulted in 
AGN with iboi/^Edd < 0.002. This accretion rate is far too 
low to fuel the obscured AGN at z ~ 1 discovered by the deep 
X-ray surveys. 

We conclude that these relatively simple models suggest 
that in order for starburst disks to be a common source of 
AGN obscuration and fueling, an efficient angular momen- 
tum transport mechanism must operate, and a global structure, 
such as a bar or spiral arms, is most likely needed to provide 
such a mechanism. 

4.4. Starburst Disks and Unobscured AGN 

The focus of this paper has been on the properties of model 
starburst disks that may potentially both obscure and feed an 
AGN. Only 41% of the model disks satisfied the selection cri- 
terion (see Sect. I2.21 i. Here, we discuss the properties of the 
remaining 59% of models that did not meet the imposed con- 
ditions. By construction, these 748 starbursts disks do not 
have a region where the central temperature surpasses the dust 
sublimation temperature and produces an inner parsec-scale 
starburst. Therefore, the photosphere of these disks is always 
thin and will not obscure any central AGN. 

The vast majority of this population of starburst disks are 
unable to fuel a central black hole at rates necessary to pro- 
duce a luminous AGN. In fact, 54% of these disks would only 
produce an AGN with Lboi/^Edd < and 75% with 

-^boi/I'Edd < 10~ 2 . These models are ones where r* < r a d v 
and so nearly all the gas in the disk is converted to stars at dis- 
tances ~ i? ut- The SFRs are within an order of magnitude to 
those plotted in Fig. |4] but are in the cold, outer regions of the 
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disk. Observationally, such systems would appear as nuclear 
starbursts within an inactive galaxy. 

Interestingly, 25% of these 'unsuccessful' starburst disks, 
or 188 overall, can fuel an AGN sufficient to result in 
-kboi/^Edd > 10~ 2 . These are models where the central tem- 
perature rises toward small radii, but just does not obtain the 
sublimation temperature. There is no pc-scale burst of star- 
formation, but star-formation does occur at modest rates over 
the entire extent of the disk. Observationally, these models 
would manifest themselves as unobscured, or Type 1, AGN. 
However, this class of starburst disks exists in a very small 
region of the explored parameter space, indicating that very 
precise conditions are required in order to result in t* < r a d v - 
This implies that if a starburst disk is a viable mechanism to 
fuel the central black hole, an obscured AGN in more likely to 
be produced than an unobscured one. This conclusion agrees 
with the observed ratios of obscured to unobscured Seyfert - 
like AGN at z ~ 1 (e.g. JBarger eTMxM. iGilli et alJl2007h . 

Collecting together the results of all the starburst disk mod- 
els reveals a possible AGN evolutionary scenario. The results 
of Sect. 13. II show that obscured accretion will most likely oc- 
cur at large values of / gas . For small values of / gas , the most 
likely outcome of a starburst disk will be to turn the gas into 
stars at large radii and an inactive central black hole. Thus, 
galactic nuclei may spend a large fraction of their history with 
small amounts of ongoing nuclear star-formation and a dor- 
mant black hole. However, processes at larger radii associ- 
ated with the assembly of the host galaxy will stochastic send 
gas-rich material toward the nucleus. This material would 
have a large / gas and, upon passing through the starburst disk, 
would result in an obscured, luminous AGN. The gas supply 
would slowly be eroded, dropping / gas and removing the in- 
ner parsec-scale starburst, revealing an unobscured AGN. Fi- 
nally, the gas fraction would drop to low enough levels to cut 
off the fueling of the black hole, and the system would return 
to producing stars at large radii until the remaining gas is ex- 
hausted. While speculative, this scenario may be a common 
occurrence during a stochastic phase of galaxy assembly. 

5. CONCLUSIONS 

In this paper we proposed that nuclear starburst disks may 
be an important contributor to AGN obscuration, especially 
at z ~ 1 where there is a large population of obscured 
AGN. The analytic m odel of starburst disks developed by 
Thom pson et alJ (120051) was used to explore the properties of 
such disks and determine if they were a viable and general 
method to obscure accreting black holes. We found that a 



range of conditions can produce potentially AGN obscuring 
bursts of star-formation on pc scales. These include a high 
dust-to-gas ratio and a relatively small size scale of the disk, 
both of which seem to be consistent with current constraints. 
The pc-scale bursts of star-formation are different from tradi- 
tional star-forming regions since they are produced in a region 
with a temperature ~ 1000 K. This region is not heated by an 
AGN, but is a result of the assumption of a Q = 1 stable 
star-forming disk in a black hole potential. Thus, the observa- 
tional signature of these nuclear star-forming regions will be 
additional warm/hot dust emission above what is required by 
AGN heating. This effect may be an explanation for the ob- 
served L2-10 ]ssw/ v L v {<a /Urn) ratios seen from both local and 
high-redshift AGN. In addition, the competition between star- 
formation and gas accretion results in a natural explanation for 
the decrease in the fraction of obscured AGN with luminos- 
ity. As an AGN-obscuring nuclear starburst will consume gas 
that may have been destined for accretion onto a black hole, 
this model predicts that nuclear starbursts should be closely 
associated with lower luminosity AGN. In contrast, high lu- 
minosity AGN require a significant gas supply, so a nuclear 
starburst would be an hindrance to fueling the black hole. 

The observational signatures of these disks may be most 
easily found in future deep radio and far-IR surveys. As- 
suming the radio-IR correlations hold for the circumnuclear 
starbursts, the expected 1 .4 GHz fluxes from these disks are 
~ 10 /iJy at the redshifts where the fraction of obscured AGN 
is near its peak. The 100/im fluxes are predicted to be ~ 1— 
10 mJy, which will be detectable in future Herschel surveys. 

Nuclear starbursts may prove to be a compelling method of 
studying the symbiotic relationship between black holes and 
their host galaxies in the era of rapid evolution of both pop- 
ulations. While qualitatively interesting, this model needs to 
be explored further to make robust quantitative predictions. 
A proper estimation of the vertical structure of the disk is re- 
quired to estimate covering fractions. Future work should also 
include incorporating the feedback from the central AGN. An- 
other important future question is to determine the nature and 
observable properties of the starburst remnants. 
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